Snowpack shifts cyanobacterial community in 
biological soil crusts 
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Abstract: Winter snowpack is an important source of moisture that influences the development of 
biological soil crusts (BSCs) in desert ecosystems. Cyanobactetia are important photosynthetic organisms 
in BSCs. However, the responses of the cyanobacterial community in BSCs to snowpack, snow depth and 
melting snow ate still unknown. In this study, we investigated the cyanobacterial community composition 
and diversity in BSCs under different snow treatments (doubled snow, ambient snow and removed snow) 
and three snow stages (stage 1, snowpack; stage 2, melting snow; and stage 3, melted snow) in the 
Gurbantungeut Desert in China. In stages 1 and 2, Cyanobacteria were the dominant phylum in the bacterial 
community in the removed snow treatment, whereas Proteobacteria and Bacteroidetes were abundant in 
the bacterial communities in the ambient snow and doubled snow treatments. The relative abundances of 
Proteobacteria and Bacteroidetes increased with increasing snow depth. The relative abundances of 
Cyanobactetia and other bacterial taxa were affected mainly by soil temperature and irradiance. In stages 2 
and 3, the relative abundance of Cyanobacteria increased quickly due to the suitable soil moisture and 
irradiance conditions. Oscillatoriales, Chroococcales, Nostocales, Synechococcales and unclassified 
Cyanobacteria were detected in all the snow treatments, and the most dominant taxa were Oscillatoriales 
and Chroococcales. Various cyanobacterial taxa showed different responses to snowpack. Soil moisture and 
irradiance were the two ctitical factors shaping the cyanobacterial community structure. The snowpack 
depth and duration altered the soil surface irradiance, soil moisture and other soil properties, which 
consequently were selected for different cyanobacterial communities. Thus, local microenvironmental 
filtering (niche selection) caused by snow conditions may be a dominant process driving shifts in the 
cyanobacterial community in BSCs. 
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1 Introduction 


Winter snowfall plays a critical role in the water and energy balance in arid and semi-arid 
ecosystems (Zhao et al., 2016; Li et al., 2020), which are significantly affected by climate change. 
Small changes in temperature and precipitation may result in large changes in the amount and 
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duration of snow cover (Brooks et al., 2011). It is predicted that the snowfall amount will decrease 
under global warming, leading to thinner snowpack in winter (IPCC, 2013; Kapnick and Delworth, 
2013). However, a recent study showed that snow cover increased in winter in some regions of the 
Northern Hemisphere, especially after 2000 (Li et al., 2019). Biological soil crusts (BSCs) are the 
dominant unit of functional vegetation in arid and semi-arid regions (Rajeev et al., 2013), and they 
play critical roles in maintaining soil stability, improving soil structure, accumulating soil nutrients, 
retaining soil moisture, and reducing salt stress (Garcia-Pichel and Belnap, 1996; Hu et al., 2003a; 
Bowker and Belnap, 2004). However, BSCs are sensitive to changes in precipitation. Winter 
snowfall is an important source of moisture and influences the growth and development of BSCs 
in temperate desert areas, especially in Northwest China (Zhao et al., 2018). 

In winter and early spring, snow cover and snow-melting processes have profound effects on 
microbial community composition, microbial activity, soil respiration, and carbon fixation and 
mineralization in BSCs in desert areas (Aanderud et al., 2013; Li et al., 2019; Ren et al., 2020). 
Cyanobacteria are the primary colonizers and important components of soil photosynthetic 
communities in BSCs (Büdel et al., 2016; Wang et al., 2020), which affect nitrogen fixation, 
moisture retention, soil stabilization and organic carbon accumulation in desert soils (Redfield et 
al., 2002; Pushkareva et al., 2018). Although many researchers have investigated the photosynthetic 
physiology and adaptation of BSC organisms, the responses of cyanobacterial community 
composition and diversity in BSCs to snowpack and snow melting are not fully understood. 

Precipitation and temperature can significantly affect photoautotrophic communities in BSCs 
(Belnap and Lange, 2003). Generally, the response of the net photosynthesis in BSCs to the 
photosynthetically active photoflux density (PFD) follows a typical saturation curve. The light level 
required to saturate the PFD of BSCs is almost at or above 700 umol photos/(m*s). For 
Cyanobacteria in BSCs, a decrease in precipitation can lead to the functional disconnection of light- 
harvesting complexes from the photosystem II reaction center and can reduce cyanobacterial 
biomass or diversity (Lan et al., 2014; Fernandes et al., 2018). Some evidence indicates that the 
amount of polysaccharides in Microcoleus vaginatus increases with increasing light intensity (Ge 
et al., 2014). The net photosynthesis of BSC organisms shows a strong link to temperature. In an 
Antarctic dry valley, field and laboratory experiments demonstrated that the rates of net carbon (C) 
fixation and dark respiration in the Nostoc commune population depend strongly on irradiance and 
temperature to form N. commune mats with water contents higher than 3096 saturation (Novis et 
al., 2007). Bacterial communities in undeveloped BSCs shifted significantly from their initial 
stages, containing few Cyanobacteria, to their late stage, dominated (2796) by Nostoc, Microcoleus 
and Leptolyngbya phylotypes during a 77-d incubation with daily freeze-thaw cycles (Schmidt and 
Vimercati, 2019). 

In this study, we hypothesized that changes in snow conditions and stages (snowpack, melting 
snow and melted snow) promote soil niche differentiation (e.g., soil moisture, irradiance and soil 
temperature), which drives variations in cyanobacterial communities. To test this hypothesis, we 
investigated (1) cyanobacterial diversity in BSCs under different snow treatments and snow stages 
and (2) the important environmental factors driving the variation in cyanobacterial community 
structure in the Gurbantunggut Desert of Northwest China. 


2 Materials and methods 


2.4 Study area 


The field experiments were conducted in the central Gurbantunggut Desert (45?03'36"N, 
87*36'36"E), Xinjiang Uygur Autonomous Region, Northwest China. The Gurbantunggut Desert 
is the largest fixed and semi-fixed desert in China, and many areas of the desert are covered by 
lichen-dominated BSCs. This region has a continental, arid, temperate climate with hot and dry 
summers and cold winters (Su et al., 2016). The annual precipitation is 70-160 mm, falling 
predominantly from April to August, and the annual average potential evaporation is 2000—2800 
mm. Although precipitation is low, there is a stable snow stage in winter and early spring in the 
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Gurbantunggut Desert. Snow meltwater in early spring provides adequate moisture for BSC 
development in this area (Zhang et al., 2007). The annual average temperature is 5.0?C—5.7?C, and 
the maximum temperature is over 40.0?C. With a coverage rate of less than 3096, the natural 
vegetation is dominated by the small tree species Haloxylon ammodendron and H. persicum. The 
elevation ranges from 510 to 550 m a.s.l. (Zhang et al., 2007; Qian et al., 2008). 


2.2 Experimental design 


In October 2014, five fixed sample sites (10.0 mx10.0 m) were randomly established in the 
interdune areas among four sand dunes. Each sample site included nine plots (area of 1.5 mx1.5 m 
for each plot). Triplicate plots were randomly established for each of three snow treatments, 
including removed snow, ambient snow (natural condition) and doubled snow (Fig. 1). A 1-m buffer 
belt was maintained between plots to prevent the different snow treatments from influencing each 
other. The plots were established to avoid the inclusion of shrubs. For the plots in the removed 
snow treatment, transparent pieces of polycarbonate (PC) (2.0 mx2.0 m) were installed 0.6 m above 
the plot surface. The light transmittance of the PC was more than 9596 (according to measurements 
of irradiance by LI-190 photoquantum sensors (Li-COR, USA)). After each snowfall, the snow on 
the PC board was removed and added to the plots of the doubled snow treatment. The PC boards 
were removed in spring, when there was no longer any snow. 
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Fig. 1 Location of the study area in the Gurbantunggut Desert (a) and demonstration of the experimental design 
(b). DS, double snow; AS, ambient snow; RS, removed snow. 


2.3 Sampling methods 


BSC samples were collected on 15 January (stage 1, snowpack), 1 March (stage 2, melting snow) 
and 15 March (stage 3, melted snow) in 2016. The snowpack stage was determined as the time at 
which a steady snowpack had formed. The snow melting stage began when the daytime temperature 
rose above 0.0°C. Stage 1 lasted from 10 December to 26 February, and stage 2 lasted from 27 
February to 13 March. Stage 3 started on 13 March. To minimize spatial heterogeneity at each site, 
we collected triplicate BSC samples with the same snow treatment at each stage from the 0-2 cm 
layer and pooled to form one composite sample. Each snow treatment included five replicate 
samples (from five sample sites) for each stage. In total, 45 BSC samples were collected from the 
three snow treatments at the three snow stages. Each BSC sample was divided into two subsamples: 
one part of the soil sample was stored at —20?C for genomic DNA extraction, and the other part 
was air-dried for the measurements of soil physicochemical properties. 
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2.4 Soil physicochemical property analysis 


Plant materials in the soil were removed by sieving with a 2-mm pore mesh. The soil organic carbon 
(SOC) was determined by the K2Cr207 method (Walkley-Black). The total nitrogen (N) was 
determined by the CuSO4-Se powder diffusion method. The total phosphorus (P) was determined 
by the NaOH fusion-Mo Te Sc colorimetry method. The available P was determined by the 0.5 
mol/L NaHCO; leaching-Mo Te Sc colorimetry method. The available N was determined by the 
alkali hydrolyzation-diffusion method. Soil pH and electrical conductivity were measured in a 1:5 
mixture of soil and deionized water using a PHS-3C digital pH meter and a DDS-307A conductivity 
meter (Precision and Scientific Corp, Shanghai, China), respectively. Total salts were determined 
by the weight method (Chen et al., 2007; Zhang et al., 2015). 


2.5 Measurements of soil moisture and photosynthetic effective radiation 


Five soil moisture sensors (Decagon Devices, Pullman, WA, USA) were installed horizontally at 2 
cm below the soil surface in each plot. They were connected to a Decagon EM 50 datalogger 
(METER, USA) that recorded the soil volumetric water content and temperature every 15 min from 
November 2015 to March 2016. To measure the photosynthetic effective radiation (irradiance), we 
installed five LI-190 photoquantum sensors (Li-COR, USA) connected to a Li-1500 datalogger 
(Shenzhen AKS Technology Co., China) on the surfaces of BSCs in each plot. 


2.6 Soil DNA (deoxyribonucleic acid) extraction and Miseq sequencing 


Soil genomic DNA was extracted from 0.5 g soil with a MOBIO Power Soil DNA Isolation Kit 
(MOBIO Laboratories, Carlsbad, CA, USA) using standard kit protocols. The universal forward 
primer 515F (5-GTGCCAGCMGCCGCGGTAA-3) and 909R reverse primer (5'- 
CCCCGYCAATTCMTTTRAGT-3') were used to amplify the V4 hypervariable region of the 16S 
rRNA gene (Xu et al., 2020). A unique 12 bp barcode was added at the 5'-end of the reverse primer, 
and each sample had its own unique barcode in Polymerase Chain Reaction (PCR). The PCR 
mixture (25 uL) contained 1xPCR buffer, 1.5 mM MgCl, 0.4 mM of each deoxynucleoside 
triphosphate, 1.0 uM of each primer, and 0.5 U Ex Taq (Takara Bio) and 10 ng soil genomic DNA. 
The PCR amplification program included initial denaturation at 94?C for 3 min, followed by 30 
cycles of 94°C for 40 s, 56°C for 60 s, and 72°C for 60 s, and a final extension at 72°C for 10 min. 
The PCR product size was approximately 394 bp. 

Two PCRs were conducted for each sample and combined after PCR amplification. PCR 
products were purified using a SanPrep DNA Gel Extraction Kit (Sangon Biotech, China). All the 
purified PCR products were mixed in equimolar amounts, used for library construction with a 
TruSeq DNA kit, and sequenced using an Illumina MiSeq sequencer. 


2.7 Sequence data analysis 


Sequence processing, clustering, taxonomic assignment and biodiversity calculations were 
performed with the QIIME pipeline (http://qiime.org/). First, sequences were demultiplexed, and 
primer and barcode sequences were removed. The sequences with high quality (length higher than 
350 bp, without ambiguous 'N' bases, and average base quality score higher than 30) were used for 
downstream analyses. To compare the relative differences among samples, we selected 15,770 
sequences per sample for bacteria randomly according to the rarefaction curve of the bacterial 
sequence data. Operational taxonomic units (OTUs) were generated by an open-reference OTU 
selection protocol, where sequences were clustered against the Greengenes database at a 97% 
similarity cutoff. Taxonomic assignments were performed using the Ribosomal Database Project 
(RDP) Classifier program (http: //rdp.cme.msu.edu/) with a confidence cutoff of 80%. The original 
sequencing reads were deposited into the National Center for Biotechnology Information (NCBI) 
database (SRA accession: PRJNA490878). 

Based on the taxonomic results, we resampled the cyanobacterial sequences from the bacterial 
sequences. To further determine the composition of cyanobacterial species, we selected 
representative sequences of the first 125 OTUs and blasted against the NCBI database. Then, 158 
reference cyanobacterial sequences of known or higher taxonomic ranks from BSCs in deserts were 
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selected. The representative sequences of the OTUs in this study and the reference cyanobacterial 
sequences were aligned using CLUSTALX (1.83) and trimmed at each end to the same length. A 
neighbor-joining tree was constructed for the aligned sequences using the maximum composite 
likelihood model in MEGA 5.0 (Auckland, New Zealand). According to the phylogenetic tree (Fig. 
S1), we chose the highest Basic Local Alignment Search Tool match to the reference cyanobacterial 
species to represent the taxonomy of the OTU (Zhang et al., 2016). 


2.8 Statistical analysis 


We calculated the cyanobacterial community dissimilarity based on Bray-Curtis distance using 
OTU tables as the input data (Rui et al., 2015). The OTU number (based on a 9796 similarity cutoff) 
and alpha diversity metrics (Observed-species, Chao 1) were produced for each sample in 
Quantitative Insights Into Microbial Ecology using the OTU table (Bokulich and Mills, 2013). The 
relative abundances of cyanobacteria at the phylum, order and OTU levels were determined by the 
percentage of the number of cyanobacterial sequences to the total bacterial sequence reads. 
Repeated-measures ANOVA (SPSS 13.0) was performed to test the differences in soil 
physicochemical properties (organic C, total N, total P, total K, available N, available P, available 
K, pH, electrical conductivity, total salts, soil moisture, soil temperature and irradiance), diversity 
indices (observed OTUs and Shannon-Wiener diversity) and relative abundances at the order and 
OTU levels among the different snow treatments and stages. When significant differences in 
different treatments and stages were observed, multiple comparisons were performed with the least 
significant difference (LSD) test. Pearson's correlation coefficients (two-tailed) among 
physicochemical properties were calculated using SPSS 13.0 prior to canonical redundancy 
analysis (RDA) to eliminate those variables that might distort the results due to their 
autocorrelations, at P<0.01 level. Considering the results of the correlation analysis, organic C and 
total K were chosen to represent the nutrient group, while pH, total salts and electrical conductivity 
were selected to represent the soil salinity group. Soil moisture, soil temperature and irradiance 
were also included as meteorological factors. 

The relative abundances of Cyanobacteria were square root-transformed prior to the analyses. 
PERMANOVA tests based on the Bray-Curtis distance measures were performed to test the 
differences in cyanobacterial community composition among the different snow treatments and 
stages using PAST (http://folk.uio.no/ohammer/past/). The correlations among cyanobacterial 
community dissimilarity and environmental variables were evaluated using the (partial) Mantel test 
with PASSaGE (http://www.passagesoftware.net/). To explore the relationships between the 
cyanobacterial community (based on OTU data) and environmental factors, we performed RDA 
with the rda function using the vegan package with the selected environmental variables and 
Hellinger-transformed data from the cyanobacterial OTU tables as the input data. The significance 
of environmental variables and RDA axes were determined by the envfit function using vegan in R 
(Borcard et al., 2011). 


3 Results 


3.1 Soil physicochemical properties under different snow treatments 


In the Gurbantunggut Desert, the ambient snow treatment received approximately 250 (+50) mm 
(n=15) of snowfall in the winter of 2015 and early spring of 2016. The ambient snow treatment also 
experienced some periods that were relatively free from snow, especially at the end of the winter 
season. Stage 1 was from 10 December 2015 to 26 February 2016. During this stage, the soil 
temperature ranged from —6.3?C to —2.9°C in the doubled snow treatment, from —7.5°C to —2.8°C in 
the ambient snow treatment, and from —22.7°C to 5.4?C in the removed snow treatment. No freeze- 
melt cycles were observed in any of the three snow treatments. In January, the snow depths in the 
ambient snow and doubled snow treatments were approximately 250 and 500 mm, respectively, 
while the removed snow treatment received little snow. Stage 2 was from 27 February 2016 to 13 
March 2016; the soil temperature under the doubled snow, ambient snow and removed snow 
treatments varied from —5.5?C to 11.9°C, —4.9?C to 8.5?C and —8.5?C to 22.1°C, and there were 
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10, 6 and 15 freeze-melt cycles, respectively. Stage 3 began on 13 March 2016, and there were no 
occurrences of snow cover or freeze-melt cycles under the doubled snow, ambient snow and 
removed snow treatments. 

The soil moisture (P<0.01), soil temperature (P<0.01) and irradiance (P<0.01) were all 
significantly influenced by the snow treatments, snow stages and their interactions. In stages 1 and 
2, the soil moisture in the doubled snow and ambient snow treatments was significantly higher than 
that in the removed snow treatment (P<0.001). The soil moisture in the doubled snow and ambient 
snow treatments in stage 2 was significantly higher than that in stages 1 and 3 (P<0.05). The soil 
temperature and irradiance increased significantly from stage 1 to stages 2 and 3 (P<0.001). In 
stage 1, the soil temperature in the doubled snow and ambient snow treatments was significantly 
higher than that in the removed snow treatment. The irradiance in stages 1 and 2 significantly 
increased with decreasing snow depth (P<0.001). The snow treatments significantly affected SOC 
(P<0.05), total N (P<0.05), available K (P<0.01) and pH (P<0.001), while the snow stages 
significantly influenced SOC (P<0.001), total N (P<0.05) and electrical conductivity (P<0.001). 
The available P was significantly affected by the interactions of snow treatments and stages (Table 
1). The SOC and total N contents in the ambient snow and removed snow treatments increased 
from stage 1 to stages 2 and 3. In stage 1, the available P content in the ambient snow treatment 
was significantly higher than that in the removed snow treatment (P<0.05). In stages 2 and 3, the 
electrical conductivity in the ambient snow treatment was higher than that in the doubled snow 
treatment and lower than that in the removed snow treatment (Table 2). These results indicated that 
both snow depth and stage influenced the soil nutrient contents (SOC, total N, available P and K), 
pH and electrical conductivity by altering the soil moisture and temperature; these changes may 
create different niches for Cyanobacteria in BSCs. 


Table1  Repeated-measures ANOVA results for the effects of snow treatments, snow stages and their interactions 
on the environmental properties of biological soil crusts (BSCs) 


Environmental — Fsnowsreatments — Fsnowstages Fsnow Environmental ^ Fnow treatments  Fsnowstage Fsnow 
property (df-2) EUER ey (df-2) DONEC 
SOC 4.489* 5.995% 3.256* pH 10.381" 2 117 0.445 
Total N 5.219; 5.266" 4.849” EC 0.947 13.675" 2.307 
Total P 0.580 0.217 0.187 Total salts 3.468 0.667 1.031 
Total K 1.041 2.261 0.525 Soil moisture 19.179" 129.255” 5.565" 
Available N 3.290 0.580 1.979 Soil temperature — 1685.724" 37.588" 409.824" 
Available P 0.783 1.288 5.078” Irradiance 3999.179% 793.259” 409.196 
Available K 6.858" 0.315 0.612 


Note: SOC, soil organic carbon; EC, electrical conductivity. * and * indicate significant correlations at the 0.01 and 0.05 levels, 
respectively. 


Pearson's correlation analysis indicated that SOC was significantly and positively correlated with 
total N (r=0.817, P<0.01), total P (r=0.357, P «0.05), available P (r=0.490, P<0.01), available K 
(r=0.661, P<0.01) and soil temperature (r=0.430, P<0.01). In addition to SOC, the total N was 
positively correlated with available P (r=0.364, P<0.05) and soil temperature (r=0.387, P«0.01), 
and negatively correlated with available K (r= —0.525, P«0.01) (Table S1). 


3.2 Cyanobacterial composition of BSCs 


In BSC samples, the number of bacterial OTUs ranged from 18,619 to 18,810. The bacterial 
communities were composed mainly of Proteobacteria and Cyanobacteria, followed by 
Actinobacteria, Bacteroidetes and Acidobacteria (Fig. S2). In stages 1 and 2, Cyanobacteria were 
predominant in the bacterial community in the removed snow treatment and less abundant in the 
doubled snow and ambient snow treatments. Additionally, the relative abundance of Proteobacteria 
in soils was significantly higher (P«0.05) in the ambient snow and doubled snow treatments than 
in the removed snow treatment. Bacteroidetes had the highest relative abundance in the doubled 
snow treatment, followed by the ambient snow treatment, and it was rare in the removed snow 
treatment. In stage 3, Cyanobacteria dominated the bacterial community in all three snow 
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treatments, and Bacteroidetes made a very small proportion of the community, suggesting the 
possible filtering role of irradiance within the microbial community (Fig. S2). 

At the order level, the cyanobacterial communities included Oscillatoriales (212 OTUs), 
Chroococcales (167 OTUs), Nostocales (32 OTUs), Synechococcales (11 OTUs) and unclassified 
Cyanobacteria (165 OTUs) in all the treatments. Oscillatoriales and Chroococcales were the 
dominant taxa, accounting for 31.1%-48.7% and 22.9%-38.3% of all cyanobacterial sequences, 
respectively. 

At the OTU level, 796 cyanobacterial OTUs were detected in BSCs. The dominant 
cyanobacterial OTUs in the different snow treatments were assigned to Microcoleaceae, 
Coleofasciculaceae, Phormidaceae, Chroococcidiopsidaceae, Chroococcaceae, Nostocaceae, 
Scytonemataceae, Tolypothrichaceae and Leptolynbyaceae. Sequences related to M. vaginatus and 
Chroococcidiopsis sp. were dominant in the BSCs, followed by M. steenstrupii, Chroococcidiopsis 
thermalis and Trichodesmium contortum (Table S2). Some OTUS related to Streptophyta were also 
found among the Cyanobacteria. A considerable proportion of the sequences could not be 
confidently classified at the genus level. 


3.3 Effects of snowpack on cyanobacterial diversity and relative abundance 


The Shannon-Wiener diversity of Cyanobacteria did not show significant changes among snow 
treatments (P>0.05) or stages (P 0.05). In stage 1, the relative abundance of Cyanobacteria was 
significantly higher in the removed snow treatment than in the ambient snow and doubled snow 
treatments (P «0.05). Additionally, the number of cyanobacterial OTUs was significantly lower in 
the doubled snow treatment than in the removed snow treatment (P «0.05; Fig. 2). When the snow 
began to melt in spring, significant increases were observed in the number of cyanobacterial OTUs 
in the doubled snow treatment and in the relative abundance of Cyanobacteria in the ambient snow 
treatment (P«0.05). When the snow had melted completely, the relative abundance of 
Cyanobacteria in the doubled snow treatment increased from stage 2 to stage 3, and no significant 
differences were observed among the three snow treatments in stage 3 (P>0.05; Fig. 2). 
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Fig.2 Effects of snow depth and stages on (a) cyanobacterial OTUs (operational taxonomic units; richness) and 
(b) relative abundance in the total bacterial community. Stage 1, snowpack; Stage 2, melting snow; Stage 3, melted 
snow; DS, double snow; AS, ambient snow; RS, removed snow. Bars mean standard deviations. The different 
lowercase letters indicate significant differences at the P«0.05 level among different snow treatments in the same 
stage (LSD). The different uppercase letters indicate significant differences at the P«0.05 level (LSD) in different 
stages of the same snow treatment. 


The repeated-measures ANOVA showed that the relative abundances of Oscillatoriales 
(P«0.05), Chroococcales (P<0.01), Synechococcales (P<0.05) and unclassified Cyanobacteria 
(P<0.05) changed significantly among the different snow treatments. The relative abundances of 
Nostocales (P<0.05) and Oscillatoriales (P «0.01) were significantly influenced by snow stage 
(Table 3). The interactions of snow treatments and stages significantly affected the relative 
abundances of Oscillatoriales (P «0.01) and unclassified Cyanobacteria (P «0.01; Table 3). In stage 
1, the relative abundances of Oscillatoriales, Chroococcales, Nostocales and unclassified 
Cyanobacteria in the removed snow treatment significantly increased (P «0.05) with decreasing 
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snow depth. In stage 2, the relative abundances of Oscillatoriales and Chroococcales increased 
significantly among the three snow treatments. In stage 3, the relative abundances of most 
cyanobacterial orders were not significantly different among snow treatments. However, the 
relative abundance of Chroococcales in the removed snow treatment was significantly higher than 
those in the doubled snow and ambient snow treatments. The relative abundances of Oscillatoriales 
in the doubled snow and ambient snow treatments increased from stage 1 to stages 2 and 3, while 
no significant differences were observed in the relative abundances of Chroococcales, Nostocales, 
Synechococcales and unclassified Cyanobacteria among the different snow stages (Table 4). 


Table3 Repeated-measures ANOVA results for the effects of snow treatments, snow stages and their interactions 
on the relative abundances of cyanobacterial orders in BSCs 


df Nostocales Oscillatoriales Chroococcales Synechococcales ius ae 
F snow treatments 2 2.364 6.193" 6.674” 4.172* 3.223" 
Fsnow stages 2 4.098" 3.798" 3.395 1.231 1.646 
Fsnow treatmentsxsnow stages 4 1.630 3.798" 2.656 0.514 2.911" 


Note: ™ and * indicate significant correlations at the 0.01 and 0.05 levels, respectively. 


In stages 1 and 2, the relative abundances of most cyanobacterial OTUs in the removed snow 
treatment were significantly higher than those in the doubled snow and ambient snow treatments 
(P<0.05) and increased with decreasing snow depth, such as those of Microcoleus vaginatus, M. 
steenstrupii, T. contortum, Symplocastrum torsivum, Symplocastrum californicum, Phormidiaceae 
cyanobacterium, Chroococcidiopsis sp. CC1, C. thermalis and Chroococcidiopsis sp. CC3 (Table 
$2). However, in stage 3, the relative abundances of all cyanobacterial OTUs did not show 
significant differences among snow treatments (Table S2). The snow stages also showed significant 
effects on the relative abundances of some species. For example, the relative abundances of M. 
vaginatus and unknown Cyanobacteria in the doubled snow treatment showed a significant increase 
from stage 1 to stages 2 and 3. In the removed snow treatment, the relative abundance of 
Coleofasciculus chthonoplastes in stage 2 was significantly higher than that in stage 1 (Table S2). 


3.4 Relationships between cyanobacterial community composition and environmental 
factors of snowpack and melting processes 


According to PERMANOVA analysis, the cyanobacterial community compositions in the doubled 
snow and removed snow treatments were significantly different (P=0.025) in stage 1. However, in 
stage 3, no significant differences were observed among the three treatments. The cyanobacterial 
community in the doubled snow treatment in stage 1 was significantly different from those in the 
removed snow treatment in stage 2 (P<0.05) and in all three treatments in stage 3 (P<0.01). 

A total of 34.86% of the cumulative variation in the cyanobacterial community was explained in 
the RDA. The first and second axes explained 18.47% and 10.03% of the cumulative variation 
(Fig. 3), respectively. The correlation analysis indicated that the first axis was significantly 
correlated with soil moisture (R?=0.261, P=0.001) and irradiance (R?-0.184, P=0.005), while the 
second axis was significantly correlated with SOC (R?=0.180, P=0.001) and pH (R?=0.114, 
P=0.016). In stage 1, obvious snow depth effects on the cyanobacterial community compositions 
of BSCs in the doubled snow, ambient snow and removed snow treatments were observed, and 
these effects were correlated with irradiance, soil moisture and pH. In stage 2, the removed snow 
treatment samples were separated from the doubled snow treatment samples. The cyanobacterial 
community compositions were correlated with soil moisture and irradiance. Most samples in the 
doubled snow and ambient snow treatments in stage 2 were also positively associated with the 
SOC and total K contents, indicating the vital role of Cyanobacteria in C fixation and nutrient 
accumulation. When the snow had just melted, the cyanobacterial community compositions in the 
ambient snow and removed snow treatments were similar, whereas they showed some differences 
from that in the doubled snow treatment. A Mantel test analysis suggested that the cyanobacterial 
communities in different snow treatments were positively associated with soil moisture 
(P=0.00060) and irradiance (P=0.00087). 
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Fig. 3 RDA plots showing the relationships between cyanobacterial OTU compositions and soil physicochemical 


properties. RDA, canonical redundancy analysis; DS, double snow; AS, ambient snow; RS, removed snow; SOC, 
soil organic carbon; EC, electrical conductivity. 


4 Discussion 


Cyanobacteria are important autotrophic microorganisms and are primary producers of organic C 
and N in BSCs. Snowpack leads to changes in environmental conditions, such as soil moisture, soil 
temperature, irradiance and other soil properties. The present study revealed how Cyanobacteria in 
BSCs responded to snowpack, and soil moisture, irradiance and temperature were observed to be 
critical factors in shifting the structure of cyanobacterial communities. 


4.1 Role of soil moisture in determining the cyanobacterial community in BSCs 


As an important source of moisture in desert ecosystems, snow directly affects the soil moisture 
content when it melts (Zhao et al., 2018). Soil moisture showed a significant decline with reduced 
snow depth in stages 1 and 2. Additionally, the soil moisture levels in the doubled snow and ambient 
snow treatments in the different stages were significantly different (Table 2). In arid and semi-arid 
regions, the maximal efficiency of CO» fixation in BSCs remains relatively constant at optimal and 
greater than optimal soil water content levels but declines rapidly at lower soil water content levels 
(Lange, 2003). For Cyanobacteria, hydration is a critical trigger of photosynthetic activity, and 
hydration thresholds may determine cyanobacterial growth (Strong et al., 2013). M. vaginatus 
exhibits a form of hydrotaxis, migrating to the surface of the crust to cause visible greening during 
infrequent rainfall events (Rajeev et al., 2013). In stage 2, the snow began to melt, and more water 
became available, which promoted cyanobacterial growth. Thus, the relative abundances of 
dominant species (e.g., M. vaginatus and Chroococcidiopsis sp. CC1) in the doubled snow and 
ambient snow treatments in stage 2 were significantly higher than those in the same treatments in 
stage 1 (Table S2). Additionally, many other cyanobacterial OTUs showed a similar pattern to the 
dominant species, although those changes did not reach a significant level. As a result, the relative 
abundances of Oscillatoriales in the doubled snow and ambient snow treatments increased 
significantly from stage 1 to stage 2 (Table 4). This could be partially explained by the remarkable 
differences in soil moisture among the different snow treatments and stages. The differences in soil 
moisture among different snow treatments and stages led to variations in the relative abundances 
of Cyanobacteria. Although desert Cyanobacteria generally have a high capacity to withstand 
desiccation (Potts, 1999; Belnap and Lange, 2003), our data suggest that various cyanobacterial 
taxa respond differently to changes in soil moisture caused by different snow depths and stages and 
that these changes can result in an increase or decrease in the relative abundances of certain 
cyanobacterial taxa under different snow conditions. 


4.2 Effects of irradiance on the cyanobacterial community in BSCs 


The relative abundances of most cyanobacterial OTUs in the removed snow treatment in stage 1 
were significantly higher than those in the ambient snow and doubled snow treatments (P«0.05). 
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In stage 2, some cyanobacterial OTUs showed similar trends as in stage 1 (Table S2). Irradiance is 
a critical factor enhancing cyanobacterial growth under extremely dry conditions (Kiderova and 
Elster, 2019). As the photic zone is generally limited to a few millimeters beneath the soil surface, 
Cyanobacteria are significantly more abundant in the BSCs than in deeper soil layers (Garcial- 
Pichel et al., 2013; Steven et al., 2013). In stage 1, irradiance increased with decreasing snow depth 
(P<0.05; Table 2). More light reached the Cyanobacteria in the removed snow treatment, which 
promoted cyanobacterial growth and reproduction. In contrast, in winter, the Cyanobacteria in the 
ambient snow and doubled snow treatments received little light (<55.8 umol/(m*-s)), and their 
growth was seriously inhibited. This resulted in lower relative abundances of Cyanobacteria in the 
doubled snow and ambient snow treatments than in the removed snow treatment in stages 1 and 2. 
In other ecosystems, such as the Tibetan Plateau and wetlands, Proteobacteria and Bacteroidetes 
are dominant members of the bacterial community in both the upper and deeper soil layers (Guo et 
al., 2015; Ade et al., 2018; An et al., 2019). The relative abundance of Proteobacteria ranges from 
23% to 29% and is negatively correlated with soil depth (Wu et al., 2017; Yang et al., 2019). In 
alpine meadows, the relative abundance of Bacteroidetes shows a positive correlation with snow 
depth (Ade et al., 2018). In our results, Proteobacteria and Bacteroidetes in the doubled snow and 
ambient snow treatments dominated the bacterial communities in stages 1 and 2, but their 
abundance decreased sharply in stage 3, when the snow melted completely. These results indicated 
that irradiance is not necessary for the growth of Proteobacteria or Bacteroidetes. 


4.3 Regulatory effect of temperature on the microbial community 


Due to the thermal insulation provided by snow in winter, soil temperature regimes are always 
more moderate than air temperature regimes (Brooks et al., 2011). Our results indicated that snow 
cover increased the soil temperature in winter. When the snow melted in spring, the soil temperature 
in the doubled snow and ambient snow treatments was significantly lower than that in the removed 
snow treatment (Table 2). Soil temperature is a key factor driving microbial biomass, abundance 
and community structure (Yuan et al., 2012; Guo et al., 2015; Vries and Griffiths, 2018). In stages 
1 and 2, the higher soil temperature and lower irradiance in the doubled snow and ambient snow 
treatments may have created unique niches to which Proteobacteria and Bacteroidetes may be better 
adapted than other taxa. Low soil temperatures and freezing conditions may lead to microbial 
dormancy or death (Ade et al., 2018). The soil temperature in the removed snow treatment in stage 
1 was very low (—11.06°C (+0.15°C)), which may have restricted the growth of bacteria such as 
Proteobacteria and Bacteroidetes. However, Cyanobacteria can withstand extremely low 
temperatures. Viable Cyanobacteria are often detected in permafrost soils and within or on the 
exterior of rocks in the desert landscapes of polar regions (Vincent, 2007). Continental-scale survey 
of the BSC cyanobacterial communities across arid and semi-arid regions indicated that M. 
vaginatus is more thermotolerant than M. steenstrupii (Garcia-Pichel et al., 2013). The ability of 
Cyanobacteria to adapt to low temperatures partially explains why the relative abundance of 
Cyanobacteria in the removed snow treatment in stage 1 was higher than those in the doubled snow 
and ambient snow treatments (Fig. 2). In stage 3, when the snow melted completely, the soil 
temperature increased, and irradiance was no longer a limiting environmental factor. Under these 
circumstances, Cyanobacteria became dominant in bacterial communities under all snow 
treatments. 

Climate warming leads to altered precipitation patterns, increased temperatures, and reduced 
snow cover (IPCC, 2013). A change in snow cover can directly or indirectly affect the soil bacterial 
community by, for example, promoting the dominance of heterotrophic Proteobacteria and 
Bacteroidetes or the dominance of autotrophic Cyanobacteria. We speculate that the long-term 
changes in snow cover induced by climate change may alter the dynamics of various microbial 
populations and their associated functions and may ultimately alter the dynamics of C and N 
processes by, for example, increasing heterotrophic respiration and reducing C sequestration. 
Further research is necessary to understand the effects of long-term changes in snow cover on the 
microbial community and its associated functions in BSCs. 
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5 Conclusions 


In summary, this study revealed that changes in snowpack significantly altered bacterial community 
structures in BSCs of desert ecosystems. Cyanobacteria were dominant in the bacterial community 
of BSCs under snow removal conditions, whereas Proteobacteria and Bacteroidetes were relatively 
abundant under snow cover. Soil moisture, temperature and irradiance were critical factors shaping 
the cyanobacterial community structure in the BSCs. Cyanobacterial abundance increased quickly 
owing to the adequate moisture and irradiance conditions in the snow melting and melted stages. 
The variations in snow depth in the different snow stages promoted the formation of new niches by 
influencing the soil physicochemical properties (e.g., soil moisture, nutrients, temperature, pH, etc.) 
and irradiance, which led to changes in the cyanobacterial community of BSCs. Thus, local 
microenvironmental filtering (niche selection) caused by snow conditions may constitute a 
dominant process driving the variations in cyanobacterial communities in BSCs. 
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Fig. S1 Phylogenetic trees based on the 16S rRNA (ribosomal ribonucleic acid) gene representative sequences of 
cyanobacterial operational taxonomic units (OTUs) from biological soil crusts (BSCs) samples and reference 
sequences using maximum likelihood method 
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Fig. S2 Bacterial community compositions in different snow treatments and snow stages. Stage 1, snowpact; 
Stage 2, melting snow; Stage 3, melted snow; DS, double snow; AS, ambient snow; RS, removed snow. The data 
in the figure are the mean relative abundance of five replicate samples. The data at top of each bar is the total 
numbers of bacterial OTUs. 
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